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ZERO TEMPERATURE LIMIT FOR INTERACTING BROWNIAN 
PARTICLES. I. MOTION OF A SINGLE BODY 1 

By Tadahisa Funaki 
University of Tokyo 

We consider a system of interacting Brownian particles in R d 
with a pairwise potential, which is radially symmetric, of finite range 
and attains a unique minimum when the distance of two particles be- 
comes a > 0. The asymptotic behavior of the system is studied under 
the zero temperature limit from both microscopic and macroscopic 
aspects. If the system is rigidly crystallized, namely if the particles 
are rigidly arranged in an equal distance a, the crystallization is kept 
under the evolution in macroscopic time scale. Then, assuming that 
the crystal has a definite limit shape under a macroscopic spatial 
scaling, the translational and rotational motions of such shape are 
characterized. 

1. Introduction. This paper is concerned with a certain scaling limit for 
a finite, but very large, system of interacting Brownian particles in M. d . The 
positions of N particles at time t, which are denoted by x(t) = (xi(t))^ =1 G 
(M. d ) N , evolve according to the stochastic differential equation (SDE) 

(1.1) dx i (t) = ~V Xi H(x.(t))dt + dw i (t) 1 l<i<N, 

where (wi(t))^L 1 is a family of independent d-dimensional standard Brown- 
ian motions. The Hamiltonian H (x) of the configuration x = {xi)f =l G (IR 6 ')^ 
is defined as a sum of pairwise interactions between particles: 

(1.2) ff(x)= J2 U{xi-Xj). 

l<i<j<N 

The potential U = U(x),x G M. d , is radially symmetric, smooth, of finite 
range and has a unique nondegenerate minimum at |x| = a > 0; see Assump- 
tion I stated in Section 2 for details. The gradient V Xi H(x) = J2j^i (xj — 
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Xj ) G R is taken in the variable . The parameter (3 > represents the in- 
verse temperature of the system. 

The basic scaling parameter e > is the ratio of the microscopic spa- 
tial unit length to the macroscopic one. The configuration x = (xi)^L 1 is a 
microscopic object and its macroscopic correspondence is given by (exi)f =l 
under the spatial scaling x i— ► ex. The goal of this paper is to investigate 
the asymptotic behavior as e j of the system defined by (1.1) with N, 
(3 and t suitably scaled depending on e, especially when the temperature 
of the system converges sufficiently fast to 0. The time change from the 
microscopic to the macroscopic levels will be introduced for x(i) by 

(1.3) x^(t) =x(e" K t), t>0, K = d + 2. 

We say a rigid crystal is formed at the microscopic level, if the particles are 
arranged in an equal distance a and the total energy H (x) increases under 
any deformation for such arrangement except isometric transformations; see 
Section 2. As j 0, that is, under the zero temperature limit, the system 
is expected to be frozen and rigidly crystallized. 

The results of this paper are twofold and will be formulated at both the 
microscopic and macroscopic levels in space. The result at the microscopic 
level can be roughly stated as follows. If the initial configuration x(0) of 
the system is nearly a rigid crystal, so is for x( £ )(t) asymptotically with 
probability one as e { if the temperature of the system decreases to 
sufficiently fast; see Theorem 3.4. 

The motion of the crystal at the macroscopic level is observed under the 
spatial scaling x \— > ex. Assuming that the particles' number N = N(e) be- 
haves as pe~ d with a fixed p > and the crystal has a limit density p(y),y G 
R as e i under the spatial scaling at time t = 0, we shall prove that x( £ )(£) 
also has a limit density pt(y) for t > 0, which actually coincides with the ini- 
tial density being isometrically transformed so that pt(y) = Pi'Pg^ v ( t ) (2/)) 
for some 6{t) and i](t). Here, (fg „ denotes an isometry on R d defined by 
<PoM = Oy + rj, y G R d for 9 = (9«%p =1 G SO(d), V = ( V a ) d a=1 G R d ; SO(d) 
stands for the d-dimensional special orthogonal group. In other words, the 
macroscopic limit ofx( £ )(£) is a rigid body with density p(y), which is congru- 
ent to the initial body. The translational and rotational motions (j](t),9(t)) 
of the limit body are random and mutually independent. They are charac- 
terized as follows: 

rj(t) = (d-dimensional Brownian motion)/-^, 

while 9(t) is a Brownian motion on SO(d) which is a solution of an SDE of 
Stratonovich's type 



d9(t) = 9(t)odm(t). 
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Here m(t) = {m a ^{t)) < ^ (3=1 satisfies m a ^(t) = —m l3a (t) and the upper half 
components {m a ^ (t);a < /?} of the matrix m(t) are mutually independent 
such that 

m a P(t) = (one-dimensional Brownian motion) /Vq + (f 

with q 01 = J Rd (y a ) 2 p{y) dy, when the coordinate y = (y a )a=i of M rf is cho- 
sen in such a manner that j Rd yp(y) dy = and (J Rd y a y^p(y) dy)^ (9=1 is a 
diagonal matrix; see Theorem 4.3 and Corollary 4.4. The constants p and 
q a + represent the total mass and moments of inertia [13] of the rigid 
body with density p(y), respectively; note that p = J Rd p{y)dy holds. 

In Section 2, the notion of rigid and infinitesimally rigid crystals is in- 
troduced together with several examples. Main results are formulated and 
proved in Sections 3 and 4. The reason that the time scaling (1.3) is the 
right one is easily observed for the (macroscopic) translational motion [see 
the identity (4.15) in the proof of Theorem 4.3], although it may not be 
obvious for rotation. At the microscopic level, the crystal translates much 
faster than it rotates. Section 5 is devoted to the proof of technical estimates 
which are needed for the proof of Theorem 4.3. Section 6 contains concluding 
remarks. 

One of the motivations of this paper comes from the theory of interfaces 
which appears under the phase transitions. The macroscopic body we have 
introduced can be regarded as a kind of Wulff shape (see, e.g., [3]) at temper- 
ature zero. The static theory for the Wulff shape is recently well developed. 
This paper attempts to analyze the motion of the Wulff shape by proposing 
a simple model. Indeed, at temperature zero, at least two pure phases arise 
in our model. One is the high density region where particles are arranged in 
an equal distance a and the other is the empty region where there are no 
particles. In this respect, the body with density p(y) is a mixture of high den- 
sity and empty regions observed macroscopically. It would be more natural 
and desirable to study the model with temperature being sufficiently small 
but fixed under the scaling. However, this problem turns out to be quite 
hard. Actually, to solve such a problem, we need to have information on 
the phase transition for the Gibbs measures corresponding to our dynamics 
with infinitely many particles, but it is not well known. Lang [14] considered 
a system of ordinary differential equation (1.1) dropping Brownian motions 
with (5 = 1 and N = oo. Such a system is obtained from the SDE in the zero 
temperature limit f3 — > oo under a time change t \—> (3~ 1 t. In one dimension 
and for strictly convex potential U having a hard core, ergodic properties of 
the dynamics were studied and equilibrium states (called rigid states) were 
characterized. Related problems were discussed for the stochastic partial 
differential equations in [4] and [5]. 

This paper deals with the motion of a single body (or single crystal in 
microscopic aspect). The coagulations of several bodies are discussed in [6] 
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in one dimension. The study of coagulations in higher dimensions is out of 
reach at present. 

2. Rigid configurations. We introduce the notion of rigidity and infinites- 
imal rigidity for configurations of particles in W d and expose several exam- 
ples of infinitesimally rigid configurations. The number N of particles is fixed 
throughout this section. 

2.1. Hamiltonian and rigidity of configurations. The Hamiltonian H(x) 
of x £ (M. d ) N is introduced by the formula (1.2). The potential U is assumed 
to be radially symmetric; that is, U(x) = U(\x\), for U = U(r),r > 0, 

and the function U(r) satisfies the following conditions: 

Assumption I. (i) (Smoothness, finite range). U £ Cq(R), where U(—r) := 
U(r). 

(ii) There exists a unique a > such that U(a) = min r >oC/(r) and c := 
U"{a) > 0. 

Condition (ii) means that the energy for two particles takes minimal value 
when the distance between these particles becomes a. The range of U is 
defined by b := inf{r > 0; U(s) = for every s > r}. Let z = (2^1 € (R d ) N 
be a configuration satisfying 

(2.1) \zi — Zj\ = a or \z% — Zj\>b 

for every 1 < i ^ j < N. An additional condition on b is necessary for such z 
to exist; for example, see condition (2.6) for configurations on a triangular 
lattice. The configuration z is a critical point of the Hamiltonian H. This 
physically means that z is a microscopically crystallized or frozen configu- 
ration of atoms at temperature zero. Its rotated and shifted configuration 
(pg :V (z) := (ipg tV (zi))fL 1 = {0zi + n)f =l is obviously a critical point of H again 
for every 6 £ SO(d) and r\ £ M. d . We shall write 

(2.2) M = We, v (z);8 £ SO(d),r) £ R d } C (R d ) N , 
and its tubular neighborhood 

(2.3) M 2 {5) = {x£ (M d ) JV ;dist(x,7W) <5}, 5>0, 

where the distance is defined under the Euclidean norm in (M. d ) N : dist(x, M) = 
inf yg _A4 ||x — y H2 and ||x — y || 2 = (J2iL\ \%i — Uil 2 ) 1 ^ 2 - The configuration z sat- 
isfying (2.1) will be called a crystal. 

We say the crystal z is rigid if H(x) > H(z) holds for every x £ M2(6)\M 
and some 5 > 0. The rigidity means that z has no internal degree of freedom 
except for the isometry. For example, in two dimension, the three vertices of 
an equilateral triangle form a rigid crystal, but the four vertices of a square 
do not. The rigid crystal is a local minimum of H by definition, but not 
necessarily a global one. 
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2.2. Orthogonal decomposition ofx£ M2($) • In order to study the rigid- 
ity of a crystal z, we introduce a decomposition of x £ A^2(^)- Let TC Z = 
{Xz + h; X £ 50(d), h £ R d } C (R d ) N be the tangent space to M at z, where 
Xz + h:= (Xzi + h)f =1 and so(d) = {X £ M(d);X + *X = 0} is the Lie 
algebra of SO(d). The set M(d) stands for the family of all d x d real ma- 
trices. The orthogonal subspace to 7i z in (M. d ) N under the inner product 
(h,h') := Td=i{h h h'i) for h = (hi)f =1 and h' = (h[)^ =l is denoted by H±. 
The corresponding norm of h is ||h||2 = (Si^i l^il 2 ) 1 ^ 2 defined previously. 

For every x £ A^2(^)i we denote by z(x) :=y £ M. the minimizer of ||x — 
y 1 1 2 in y £ M- If 5 > is sufficiently small, z(x) is uniquely determined and 
x £ A42(5) admits a decomposition: 

(2.4) x = z(x)+h(x), z(x)eM, h(x)G^ (x) . 

In fact, by the definition of z(x), we have 

= 



-^-\\x-ip e ux M (z(x))\\l 



u=0 



for every X £ so(d) and h £ M. d , and this implies h(x) := x — z(x) £ 7^? x v 
2.3. Hessian of H on M. Let 

Hess H (x) = ( 92 H (x)) £ M(dN) 



\dxfdXj J l<a,f3<d,l<i,j<N 

be the Hessian of H and define a quadratic form in h = (hi)fL 1 £ (R ) by 

i,j=l a,/3=l ® x ? 

Then a direct calculation yields the following at x = z. 



(h,Hesstf(x)h)=£ £ ——ptohthl 



Lemma 2.1. 

fi(h) = £ 1>z (h) := (h,Hess#(z)h) = \ - ^ - > 0, 

where the sum (i,j) is taken over all pairs {i,j} satisfying \z{ — Zj\ = a. We 
call such pairs neighboring. 



This lemma immediately shows that the Hessian of H degenerates for 
h £ 7i z . The degeneracy for (hi = h)^L 1 comes from the invariance of H 
under the translation, while that for (hi = Xzi)f =l comes from its invariance 
under the rotation. The rigidity of z follows from the nondegeneracy of the 
Hessian of H for h £ 7i z . This leads us to introduce the following notion. 
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Definition 2.1. We call the crystal z infinitesimally rigid if the quadratic 
form Hess H (z) restricted on the subspace H z is (strictly) positive definite: 

^(h) = o^heH z . 

Since H{x) =H(z) +5i(h(x))/2 + o(||h(x)|||) as ||h(x)|| 2 -> under the 
decomposition (2.4) [see (3.1)], we easily see that the infinitesimal rigidity 
implies the rigidity. 

Remark 2.1. The study of rigidity and infinitesimal rigidity for bar and 
joint frameworks has a long history; see [1], [2] and [17]. The length of bars 
is always a in our case, but such an assumption is unnecessary in a general 
theory. According to ([1], page 281), z is called rigid in M. d if, for every 
x sufficiently close to z satisfying | for every neighboring pair 

there exists an isometry ip of M. d such that Xi = ip(zi) holds for every 
i. In [2] (page 173) z is called infinitesimally rigid in R d if T z = ker dfc(z) 
holds, where T z = H z and kerd/ G (z) = {h G (R d ) N ; £i(h) = 0}, respectively, 
in our terminology. We therefore see that the definitions of rigidity and 
infinitesimal rigidity employed by these papers coincide with ours. Three 
points {pi,P2,P3} in I^ 2 sitting on a line are rigid but not infinitesimally 
rigid in R 2 when the distances between any two points are specified; see [12], 
noting that the rigidity in that paper means the infinitesimal one. This 
example is not for a crystal, but exhibits the difference in two notions. 

The rigidity of z implies the connectedness of the set z under the neigh- 
boring relation and therefore we have, under the infinitesimal rigidity, 
the spectral gap for the quadratic form £i(h) in the following sense: 

A«(z) = inf { |^|; h G || Vh|| 2 + o} > 0, 

where 

\\Vhg = J2\h-hj\ 2 . 

{id) 

This can be rewritten as 

(2.5) A( 1 )(z)||Vh||2< < S 1 (h)<c||Vh||2, h G H z . 

Note that the second inequality is obvious. 

Remark 2.2. In one dimension, a chain z = [z\ = ai)jL 1 arranged in an 
equal distance a is a rigid crystal and A^(z) = c since £i(h) = c||Vh||2. 
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2.4. Examples of infinitesimally rigid crystals. We prepare two lemmas 
before constructing several examples of infinitesimally rigid crystals. We say 
a set C = {xi}2 = Q C n < d, is an n-dimensional cell in K d if the dimension 
of the affine hull of C is n and \xi — xj \ = a for every < i < j < n. 

Lemma 2.2. Let {ej}f =1 6e a basis of M. d and set eo = 0. If vectors 
{hi}f =0 C M. d satisfy (hi — hj,et — ej) = for every < i < j < d, then there 
exists a unique X £ so(d) such that hi = Xei + ho, 1 < i < d. In particular, 
a d- dimensional cell is infinitesimally rigid. 

PROOF. The uniqueness of X is obvious, since {ei}f =1 forms a basis 
of M. d . To show the existence of X, we may assume that (ei, . . . , ea) is an up- 
per triangular matrix. In fact, by Schmidt's orthogonalization, one can find 
P € 0(d) and upper triangular matrix (e"i, . . . , id) such that (e\, . . . , e^) = 
P(e"i, . . . , id). If the conclusion holds for (ei, . . . ,id), there exists X £ so(d) 
such that P~ l Xi = Xi h l<i<d. Taking X = PXP" 1 G so(d), the conclu- 
sion is shown also for {ei}f =1 . 

Now we assume (ei,...,e^) is an upper triangular matrix and use an 
induction in d to construct X. We may further assume ho = by replacing 
hi with hi — ho . Since ej has a form 



= (hi - hj,ei - ej) = (h\ - ti p e\ - e'j), l<i<j<d-l, 

and therefore h\ = X ' e\ holds for some X' £so(d— 1) and every 1 < i < d — 1 
by the assumption of the induction. Now, writing 



Then, we see that hd = Xed holds. We need to prove that hi = Xei holds 
also for 1 < i < d — 1. To this end, it is enough to show that — (y, e'j) = h'(. 




with e[ E R d 1 for 1 < i < d — 1 , writing 




with h'i G and h'( G R, we have 





d-l 
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However, since (hi - h d ,ei - e d ) = (K,ei) = (hd,ed) = 0, we have (hi,e d ) + 
(hd,ei) = 0. This implies — (y, e'j) = h![ , since (h d ,ei) = (h'^e^) and 

(h h e d ) = (h\^ d ) + h'le d = {X'4,e> d ) + fcJVj 

= -(e'„lV,) + Jtfeg = -(4, h' d ) + egCeJ, y) + h'^. 

Therefore, the conclusion is shown in d dimension if it is true in d — 1 
dimension. The procedure of the induction is complete once we can show 
the conclusion when d = 2. However, this is already essentially done in the 
above argument. 

The infinitesimal rigidity of d-dimensional cell C = {xi}f =0 is immediate 
by taking e, = X{ — xq for 1 < i < d. □ 

Lemma 2.3. Let two infinitesimally rigid crystals z^ and zP^ be given 
and assume that the dimension of the affine hull of z^ 1 ) n z^ is at least 
d—1. Then, the joined configuration z^ 1 ) Uz' 2 ' is infinitesimally rigid. 

PROOF. Let us denote zW n z^ = {zf^}^ v Then, the conclusion fol- 
lows if one can show that Xzf^ + h = X'zf^ + h' for every 1 < i < Nq implies 
X = X' and h = h' , where X, X' G so(d) and h, h' G M d . However, from the 
assumption, one can find at least d—1 linearly independent vectors {e/c}^z} 
from {z^ — zjj }i<i<j<N - The identities Xe k = X'ek hold for such vectors. 
Take e d G M d in such a manner that {e k }f =1 forms a basis of M. d . Then, since 

(Xe d ,e k ) = -(e d ,Xe k ) = -(e d ,X'e k ) = (X'e d ,e k ) 

for 1 < k < d — 1 and (Xe d ,e d ) = = (X'e d ,e d ), we see Xe d = X'e d . We 
accordingly have Xe k = X'e k for every 1 < k < d. This proves X = X' and 
therefore h = h! . □ 

Example 2.1. The set obtained by patching together <i-dimensional 
cells on their faces is infinitesimally rigid from Lemmas 2.2 and 2.3. More 
precisely, a finite set icM (i satisfying the following two conditions is in- 
finitesimally rigid: 

1. A = \J k C k with finitely many d-dimensional cells C k . 

2. For any C kl and C kn in A, there exists a sequence C k2 , . . . , C kn _ 1 in A 
such that C ki n C ki+1 are (d — 1) -dimensional cells for 1 < i < n — 1. 

In two dimension, the set as in Figure 1 is infinitesimally rigid. In general, 
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Fig. 1. Two-dimensional crystal. 



infinitesimally rigid crystals may have defects. 

Example 2.2. In three dimension, the tetrahedron ( = three-dimensional 
cell), octahedron and icosahedron are infinitesimally rigid by Cauchy's rigid- 
ity theorem or by Alexandrov's rigidity theorem [2]. Note that the faces of 
the three types of regular polyhedrons listed above are all equilateral tri- 
angles. In particular, the set obtained by patching together polyhedrons of 
these three types as in Example 2.1 is infinitesimally rigid. 

Example 2.3 (Crystals on triangular lattice). Let {e a £ R d } d =1 be a 
basis of Hl d such that (e a ,ep) = (1 + 6 a p)/2. In other words, it is a system of 
unit vectors and arbitrarily chosen two of them are at an angle of 60° with 
each other. Then, a d-dimensional triangular lattice is defined as an integer 
lattice generated by {e a } d a=1 :A = A d = {£^=1 &,e a £l d ^ = (&*)a=l e % d }. 
Note that Ai = Z. Set E = {e E A; |e| = 1}, and then it is easy to see that 
e = Xla=i £« e » S E if and only if £ has the form £ a = ±5 aQa for some ao 
or £ Q = 5 aoa — 5p a for some ao / 0o- The triangular lattice A is the set 
of centers of circles (when d = 2) or balls (when d = 3) with radius 1 /2 
filled most densely in the space; this assertion was known as the Kepler 
conjecture in three dimension and solved by Hales [8]. We need an additional 
assumption: 

(2.6) b<c(A d )a, 

for a rigid z to exist on aA satisfying (2.1), where c(A^) := ini X £A d \E \ x \ = 2 
(when d= 1), y/3 (when d = 2) and y/2 (when d > 3). 
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In two dimension, A can be constructed by patching equilateral triangles, 
while in three dimension, A is obtainable by patching tetrahedrons and 
octahedrons. Therefore, at least in two and three dimensions, e~ l D n aA 
is infinitesimally rigid for a bounded domain D in M. d (d = 2 or 3) having 
smooth boundary and for small e by deleting or adding some points near 
the boundary in a proper way if necessary. 

2.5. Tubular neighborhoods of A4 defined in two other norms. Let us 
consider two norms llhlloo 

and ||Vh||oo for h = (/tj)ili G (R d ) N defined by 
llhlloo = max \hi\ and || Vh^ = max \hi — hj\, 

1 {hi) 

respectively; recall that refers to neighboring pairs. Then, for every 

small c> 0, tubular neighborhoods Moo(c) and A^^(c) of M. can be intro- 
duced as 

(2.7) Moo(c) = {x£ (R d ) N ; ^(x)^ < c}, 

(2.8) Ml(c) = {xe (RT; l|Vh(x)|U < c}, 

respectively, where h(x) € (M. d ) N is defined by (2.4) for x € M.2($) with 
sufficiently small 5 > 0. For each crystal z, since two norms ||h||oo arm l|h||2 
are mutually equivalent, one can find c = c(z) > such that h(x) is well 
defined for all x £ A^oo(c(z)). 

For 1 < i / j < A/", let p(z, j) = {z = ~ h ~ • ■ • ~ *n = j} be the shortest 
path connecting i and j, where ifc ~ ik+i means that the pair (ik,ik+i) is 
neighboring. We call n =: (Jp(i, j) the length of p(i, j) and define the radius 
of z by 

J R(z) = max{ttp(i,j); l<i#j<iV}. 

Lemma 2.4. For every x 6 A^oo(c(z)), we Ziaue 

(2.9) ||h(x)|| 00 < J R(z)||Vh(x)|| 00 . 

In particular, the set A^^(c) is well defined for < c < c(z)/i?(z). 

PROOF. Since £)JLi = for h = h(x) = (hi)f =1 , 



1 N 



1 * 



1 N 

<]vE E |/^-/^ +1 |<^( z )ll vh (x)||oo 

i =1 («fe,«fc+i)ep(j,i) 

for every 1 < i < N. This shows (2.9) and M^(c) C Moo(i2(z)c). □ 
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3. Microscopic shape theorem. This section establishes the asymptotic 
behavior as e j of x^ e \t), the solution of the SDE (1.1) which is scaled 
macroscopically in time, when the initial configuration x( £ )(0) is nearly in- 
finitesimally rigid and the temperature j3~ l = /3(e) -1 of the system decreases 
to sufficiently fast compared with e. 

3.1. Behavior of H near M. Let z be a crystal, that is, a configuration 
satisfying the condition (2.1). A configuration x = {xi)f =l is then decom- 
posed as x = z + h around z just by setting h = (hi)f = i := (xi — zi)^L 1 . We 
shall write 

N 

G(x)=^|V^(x)| 2 (^||V J F/(x)||^). 
i=l 

Lemma 3.1. Suppose that x satisfies \xi — xj\ > b for nonneighboring 
pairs {i,j} (of z) and \xi — Xj\ > ao for neighboring pairs (i,j) with some 
ao G (0, a). Then, there exists C > (independent o/x,z and N) such that 

(3.1) \H(x) - {H(z) + i«Si(h)}| < CJ2 \hi ~ hj\ 3 , 



(3.2) \G(K)-S 2 (h)\<Cj2\hi-h j 



{hi) 

3 
(i,3) 



where 



c 2 ^ 



a . 

2=1 



(/ii - - Zj)(zi 



1 N 

i=l 



Proof. These two estimates are shown by Taylor's theorem applied for 
if (x) and G(x) in the variables {hi — hj} with neighboring pairs Since 
the computations are easy, the details are omitted. □ 

Now let us assume that z is infinitesimally rigid. 

Lemma 3.2. There exist C > and X^(z) > such that 
C- 1 S 2 (h) < Sx(h) < {X^(z)}- 1 £ 2 (h), h G H±. 

Proof. The first inequality is obvious. To show the second, note that 
the quadratic form £i(h) is expressed as £i(h) = (Ah, h) with a symmetric 
matrix A £ M(dN) and 7i. z is the eigenspace of A corresponding to the 
eigenvalue 0. Since we have 

£ 2 (h) = \\Ah\\ 2 2 = (A 2 h,h), 

£ 2 (h) = holds if and only if h G 7i z and this implies the conclusion. □ 
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Remark 3.1. The constant A^ 2 ^(z) is related to the Poincare inequality; 
see Lemma 2.1 of [6] in one dimension. We may assume 0< A( 2) (z) < 1 for 
every z. 

In the following, we shall normalize the Hamiltonian H as H (z) = by 
adding a constant [i.e., by considering H — H(z) instead of H\. 

Corollary 3.3. //h = x-ze^ and HVhUoo < 5\W{z)X^\z) is 
satisfied for sufficiently small 5 > 0, we have 

C- 1 A (2) (z)F(x) < G(x) < CH(x), 

for some C > 0. 

Proof. The estimates (3.1), (2.5) and the assumption on ||Vh||oo show 
that 

\H(x) - |£Ti(h)| < CllVhlloollVhH 2 < CSS^h), 

since A^ 2 ^(z) < 1. Therefore, taking 5 sufficiently small so that C5 < 1/4, we 
have 

(3.3) f£i(h)<ff(x)<f£i(h). 
On the other hand, from (3.2), (2.5) and Lemma 3.2, 

|G(x) - £ 2 (h)\ < CllVhlUllVhll 2 < C5£ 2 (h), 

which shows 

(3.4) |f 2 (h)<G(x)<|£2(h). 

The conclusion follows from (3.3), (3.4) and Lemma 3.2. □ 

3.2. Lyapounov argument. Assume that a sequence = (z^)^L 1 , < e< 1, 
of infinitesimally rigid and centered crystals is given, where "centered" means 
J2iLi = 0. The number N = N(e) of particles in may change depend- 
ing on the scaling parameter e. Let x(i) = (xi(t))fL 1 G (M. d ) N be the solution 
of the SDE (1.1) and introduce the time changed process x^ £ \t) of x(t) 
by (1.3). The inverse temperature (3 = (3(e) changes with e and diverges to 
+oo as e I sufficiently fast; see the condition (3.8) in Theorem 3.4. 

Let M=M^ and MZ(c) = M^ [e \c) be the sets (2.2) and (2.8) de- 
termined from z( £ ) instead of z, respectively. Given c = c(e) [ 0, consider a 
sequence of stopping times a = defined by 

(3.5) a = inf {t > 0;x^(t) $ M w J e \c(e))}. 

The main result of this section can now be stated. The proof is based on the 
Lyapounov type argument. 
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Theorem 3.4. Let {c = c(e) J 0} (as e j. 0) and a sequence of (random) 
initial data {x( £ )(0)} be given and satisfy the following conditions: 

(3.6) < c(e) < 5A (1 ' e) A (2 ' e) Ac (£) , 

(3.7) i?[||Vh(x( £ )(0))||^ 1/2p = o({A( 1 ' e )} 1/2 C ( e )), 

as e 10 for some p > 1, where = \W(z&), A^ = A^(z (£) ), c& = 

c(z^)/R(z^) and 5 > is the small constant appearing in Corollary 3.3. 
We further assume the following condition on the sequence of temperatures 

0-*= Pie)' 1 10: 

(3.8) (3(e)- 1 = ({A( 1 ' £ ) C ( e ) 2 Af(e)" 1 } p/(p - 1) A( 2 ' £ )e K /( p - 1 )), 
as e 10. Then we have, for every t>0, 

limP(cr( £ ' >t) = l. 

e|0 

Theorem 3.4 asserts that asymptotically, with probability one, x^ e \t) 
keeps its rigidly crystallized shape within fluctuations c(e). In order to make 
the fluctuations smaller, we need better assumptions on initial data as re- 
quired in (3.7) and on the speed of convergence to of f3(e)~ l as in (3.8). 
This theorem characterizes the microscopic structure of the solutions of the 
SDE (1.1), which are scaled macroscopically in time. 

Remark 3.2. (i) Condition c(e) < c^ in (3.6) is necessary only for the 
set A4^ e \c(e)) to be well defined; recall Lemma 2.4. 

(ii) Condition (3.7) is always satisfied if x( £ )(0) G . 

(iii) The theorem covers the situation purely microscopic in space, that is, 
the case where the particles' number N is fixed and does not change with e. 
In this case, c = c(e) can be taken independently of e if it is sufficiently small 
and the condition (3.8) is satisfied if the temperature behaves as /3(e)" 1 = 
o(e«/(p-i)). 

(iv) The result will be reformulated in one dimension in [6], Theorem 2.2, 
and the condition (3.8) will be rewritten into much simpler form based on 
Remarks 2.2 and 3.1 on A^ 1 ^ and X^ 2 \ respectively. 

For the proof of the theorem, we first note that x( £ )(i) = (x[ e \t))fL 1 sat- 
isfies the following SDE: 

(3.9) dx\ e \t) = -^e- K V Xi H(^ £ \t))dt + e- K / 2 d Wi (t), l<i<N, 

in law's sense. Simple application of Ito's formula shows the next lemma. 
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Lemma 3.5. For every p> 1, 

H p {x^{t)) = H p {^{0))+m p £ \t) 

+ A-/fe~ K Mx (£) (*)) + e- K b 2 , p (^\s))}ds, 
Jo 

where 

6 1 , p (x) = |^- 1 (x)G(x), 

N 

6 2)P (x) = |(P- l)^" 2 (x)G(x) + |if^ 1 (x)^A Xi F(x), 

with the Laplacian A Xi in the variable Xi, and nip (t) is the martingale 
defined by 

m £ )(t) = £ - s / 2 pV / ^- 1 (x( £ )( s ))(V XiJ ff(x( £ )( S )),d U ; J ( S )). 

We have the following bounds on the drift functions b\ p and b 2 , p - 

Lemma 3.6. Assume x E .A/f^'^^A^'^At 2 '^ A c~( e )), where 5 > is the 
constant appearing in Corollary 3.3. Then, there exists C = C p > such that 

(3.10) 6i iP (x) > C- l \^ £) H p (x), 

(3.11) 6 2!P (x)<CiViJP- 1 (x). 
/n particular, 

(3.12) - /36 1)P (x) + & 2 , p (x) < C/3- p+1 iV p {A( 2 ' £ )}- p+1 . 

Proof. The lower bound (3.10) is immediate from Corollary 3.3. To 
show (3.11), note that, for every < a < c 2 , | E*Ii A x^(x)| < CN holds 
for some C > and all x E (R^)^ satisfying ci < |xj — < C2 for each (i,j). 
Then, we have (3.11) from Corollary 3.3. Finally, to show (3.12), we estimate 
choosing q= p/(p — 1) (when p^l), 

HP- 1 ^) = {Lf3- 1 N{\ i - 2 ^y l ) 1,q {L' 1 f3N~ 1 \ { - 2 ^ H p {^)) l/q 

< -(Lp- 1 N(\^r 1 ) p - 1 + -L-^N^X^HPhc), 
p q 

for every L > 0. The inequality (3.12) follows from (3.10) and (3.11) by 
choosing L > sufficiently large. □ 
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PROOF of Theorem 3.4. Lemmas 3.5 and 3.6 show 
E[H p {^ £ \t A a-))] < E[H p {^ £ \o))} + Ce- K f3- p+1 N p {X^' £ '>y p+1 t. 

However, (3.3), (2.5) and the assumption on x( £ )(0) imply that Op := 
(A^ x c(e) 2 )~ p x E[H p (^ £ \o))\ tends to as e [ 0. On the other hand, if 
x = (xi)f =l E dAd^c , (c(e)), then \hi — hj\ = c(e) for some neighboring pair 
and therefore -ff(x) > C~ 1 A( 1,£ ^c(e) 2 from (3.3) and (2.5). Accordingly 
we have 

E[H p (^ £ \t A a))} > E[H p {^ £ \a)),a <t}> (C^X^ c{e) 2 ) p P{a < t). 
Therefore, we have 

P(a <t)< C{a& + (X^ £) c{e) 2 y p e- K p- p+1 N p {X^}- p+1 t}, 
which tends to as e j 0. The constants C may change from line to line. □ 

4. Motion of a macroscopic body. In this section we shall identify the 
motion of a macroscopic body obtained in the limit under the spatial scaling 
x i— > ex as e I 0. 

4.1. Coordinate #(x). Let z = (zi)f =l be a centered infinitesimally rigid 
crystal and we fix it throughout Section 4.1. For x E M. 00 {c{z)), z(x) E M. 
is defined in Section 2.2 as the minimizer of ||x — y 1 1 2 in y G see also 
Section 2.5. Since z(x) E M., one can represent it as z(x) = ipg tTI (z) for some 
{9,i]) = (#(x),77(x)) E SO(d) x R d . The function ?y(x) defined in this way 
actually coincides with the center of the mass of x: 

1 N 

(4.1) i7(x) = -J>€R d . 

i=l 

In fact, (4.1) is seen from d||x — ipQ^iz) 1 ^/ dn a = for every 1 < a < d. On 
the other hand, the function #(x) = (^"^(x))^ (3=1 has the following property. 

Lemma 4.1. For every 1 < a, (3 < d and x £ Moo(c(z)), 

N 

(4.2) (W* (x), Vff(i)) = £ (V,^ Q/3 (x), V I( ff(x)) = 0. 

i=l 

PROOF. Step 1. For every y £ M, let ee(y) £ (R**)^ and A^(y) > 0, 
1 < I < dN, be the eigenvectors and the corresponding eigenvalues of the 
Hessian Hess-ff(y) of H at y, respectively. Recalling that y is infinites- 
imally rigid, we may assume A^(y) = for 1 < £ < d and A^(y) > for 
d + 1 < I < dN, where d := d(d + l)/2 is the dimension of the space A4 
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or 7i y . Then, the vectors (ee(y))f =l and ( e ^(y))^_j +1 span the spaces 7i y 
and TCy , respectively. Moreover, since the invariance of the Hamiltonian H: 
H((p g>ri (x)) = H(x) implies HessF(^^(y)) = (p efi H(y)ipg^ one can take 
{e^(y), A^(y)} in such a manner that 

(4.3) e e (ip e>v {y)) = ip efi (e e (y)), \e(<Pe, v {y)) = A^(y), 

for every 1 < £ < dN, (6, ?]) G 50(d) x M d and y G M 

Since x — z(x) G W^ x ^ by (2.4), setting y := z(x), x G A^oo(c(z)) can be 
decomposed as 

dN 

(4.4) x = y+ J2 CVy) 

£=J+i 

for some G R, J + 1 < I < dN. We call (y, . . . , ( dN ) £Mx R^N-d the 
Fermi coordinate of x; see [7], page 4. 

Step 2. Under the Fermi coordinate, the Hamiltonian H does not depend 
on the variable y: 

(4.5) H(x)=H(( i+ \...,( dN )- 

To see (4.5), we first note that the Fermi coordinate of (pg t71 (x) is given by 
(vWy),C J+1 ,..-,0- In fact ' by ( 4 ' 4 ) and then by ( 4 ' 3 )' 

dN \ dN 



/ dN \ uiV 

',o( 2 C%(y) J =w,T,{y) + £ C%(<^,r,(y))- 



Therefore, the invariance of implies 

H( y , ( d+1 , . . . , C d7V ) = # K„(y), • • • , C d7V ) 

under the Fermi coordinate for every (#,77) G SO(d) x R d . This shows (4.5). 

Step 3. We finally prove Vff(x) = (V^P^x))^ G ^ (x) and W Q/3 (x) G 
H z ( x ). Once these relations are established, the conclusion of the lemma is 
immediately deduced. Take an arbitrary £ G (R d ) and decompose it as £ = 



P£ + P x £, where P : (R d ) N -> ft z(x) and P x : (M^ -> ?^ (x) are orthogonal 
projections. Then 



(V^(x),0 = ^(x + <) 
an 



M=0 



(46) d 

= f- P^x + uP^) 
du 



(VP(x),P^)- 
) noting that c 

uP£, A4) = 0(u 2 ) as it ->• 0. Equation (4.6) implies Vff(x) G Hi x) . Since 



The second equality of (4.6) follows from (4.5) noting that dist(z(x) + 

,2\ 



ZERO TEMPERATURE LIMIT, I. SINGLE BODY 17 

by definition 6(x) depends only on y = z(x) :#(x) = 9(z(x)), one can simi- 
larly show V6^(x) G H z(x) . □ 

In order to identify the motion of the macroscopic body, it becomes 
necessary to calculate the derivatives of 0(x) in the variables X{. We in- 
troduce notation to give their representations. The space M(d) of d x d 
matrices is equipped with an inner product (X,Y) := Tv(X t Y) and a norm 
\X\ = \/(X, X) for X,Y £ M(d). The orthogonal projection from M (d) onto 
its subspace so(d) under this inner product is denoted by Proj, that is, 
Proj X = (X - l X) /2 for X £ M(d). For e = (e a ), e = (e a ) £ M d , determine 
the matrix e <g> e = ((e <g) e) Q/3 ) £ M(d) by (e <g) e) Q/3 = e a eP , l<a,/3<d. 
The 7-directed unit vector in R d is denoted by c 7 , 1 < 7 < d. We define 
Q(x) = (g a/3 (x)) G M(d) by 

AT AT 

(4.7) Q(x)=^2z i ®x i , that is, g Q/3 (x) = ^ zfxf , l<a,/3<d. 

t=l i=l 

Note that Q(z) is symmetric. The map {Proj o(Q(x)#(x))} _1 is the in- 
verse of Proj o(Q(x)0(x)) :50(d) 3Ih Proj(Q(x)0(x)Jf) = {Q(x)6(x)X + 
X*(Q(x) x 0(x))}/2 £ 30(d). Note that the derivative 

^(x) g r fl(x) (50(d)) = e(x){so(d)}, 



the tangent space to SO(d) at 0(x). 

TION 4.2. For 1 < 7 < c 

■(x) = 0(x){ Proj o(Q( x )0(x))}- 1 Proj {^(x)" 1 ^) 8> 



Proposition 4.2. For 1 < 7 < d and 1 < z < iV, we have 
d6 



dxj 

Proof. Since ¥?0( x ),r;(x) ( z ) i s the minimizer for the norm ||x — y||| in 
y G M. , we have 



, AT 

-^Ix.-^e^-^x)! 2 



= 



for every Y" G so(d). However, since E^Li «j = and (6(x)zj,6(y)Yzj) = 
(zj,Yzj) = for all j, this implies J^L^a^-, 6(x)Yzj) =0, which can be 
rewritten as 

(4.8) (Q(x)0(x),y) = O. 

Taking the derivative of (4.8) in and noting that dq a/3 / dx] (x) = zfS^ 1 , 
we get 

(4.9) (q( x )J^(x),y) =-(0(x)(9(x),y) = ((0(x)- 1 e 7 )®*,r) 
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for every Y <E so(d). Set X := 6>(x) l d9 /dx](x.) and then, since X G so(d), 

(4.9) shows that 

Proj (Q(x)0(x)X) = Proj {(^(x)- 1 e 7 ) 8> 
This proves the conclusion. □ 

4.2. Identification of the limit. We now discuss the limit as e j. under 
the macroscopic spatial scaling x^> y = ex for the system with the particles' 
number iV = N{e) changing with e. 

Our formulation is the following. Let SDTp = Wlp(M. d ), p > 0, be the fam- 
ily of all Radon measures p on R d satisfying p(R d ) < p. The space Tip is 
equipped with the topology determined by the weak convergence. A sequence 
x (e) _ (xj £ - ) )^ 1 , < e < 1, with N = N{e) of the system of particles in R d is 
identified under the scaling with p^ £ '(x.^) £ VJlp,p = e d N defined by 

N 

(4.10) p^(dy) = ^\^-dy) :=E d Y,5 ex i,){dy). 

i=l 

Let us assume that, as in Section 3.2, a sequence z^ = (z^ )fL t , < e < 1, 
of centered infinitesimally rigid crystals is given and satisfies the following 
three conditions: 

1. There exists R > such that \z^\ < Re" 1 for all i and e. 

2. z( £ ) has a macroscopic limit density function p(y), y € M. d , in the sense 
that p^ (z^ ; dy) =^ p(y) weakly as e j 0. 

3. (Nontriviality of the limit). The total mass of the macroscopic limit den- 
sity is positive: p := f Rd p(y) dy > 0. 

Examples of the sequence z^ will be given at the end of this section. We 
shall denote by M.^{c) the set (2.7) determined from z^ instead of z. 
The domain D := {y G M. d ;p(y) > 0} with density p(y) may be called the 
macroscopic shape of the body. The above conditions imply that D C {y G 
M. d ; \y\ < R} and lim £ | e d N(e) = p. Since z^ are centered, the body with 
density p(y) is also centered in the sense that f Rd yp{y) dy = 0. Let Q = 
(o al3 ) d t 3=1 G M{d) be the matrix defined by 

<fP= f y a yPp(y)dy. 

Then, since the matrix Q is symmetric, by rotating the body around the 
origin G M. d , we may assume that Q is diagonal with diagonal elements: 

(4.11) t = I (y a ) 2 p(y)dy, l<a<d. 
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Let x( £ \t) := x(e~ K t), n = d + 2, be the process obtained by macroscop- 
ically scaling in time the solution x(t) of the SDE (1.1) with initial con- 
figuration x(0) = z( £ ) . The spatially macroscopic scaling limit of x^ e ^ (t) is 
characterized by the following theorem, in which the limit of (n^ (t) , 9^ (t)) := 
(en(^ e \t)),9(^ £ \t))) G R d x SO(d) as e j is obtained. Here, the coordi- 
nate 0(x), x G M ( ^(c(z^)), is defined as in Section 4.1 based on z( e ) in 
place of z. If x( £ )(t) goes outside of .M^(c(z( £ ))) at a certain time, 9^ e \t) 
may be defined arbitrarily after such time keeping it continuous in t. This 
theorem, in particular, shows that the proper macroscopic time scalings for 
the translational and rotational motions of the body are the same. 

Theorem 4.3. Assume that the temperature (3~ l = /3(e) -1 of the system 
converges to as in (3.8) for some c(e) satisfying (3.6) and c(e) < e v for 
some v > 0. Then, the process (r/ £ ) (t) , 9^ £ \t)) weakly converges to (rj(t),9(t)) 
as e I in the space C([0,T],M. d x SO(d)) for every T > 0. The limit is 
characterized by the following three properties: 

(i) n(t) and 9(t) are mutually independent. 

(ii) y/pn{t) is a d-dimensional Brownian motion starting at 0. 
(hi) 9(t) is a solution of an SDE of Stratonovich' 's type on SO(d): 

(4.12) d9(t) = 9(t)odm(t), 9(0) = I, 

where m(t) = (m a ^(t)) d i g =l is an so(d)-valued Brownian motion such that 
the components {m a ^ (t);a < {3} in the upper half of the matrix m(t) are mu- 
tually independent and ^Jq\ is one- dimensional Brownian mo- 
tion for each 1 < a, (3 < d. 

This theorem can be reformulated as the convergence for measure-valued 
processes. 

Corollary 4.4. Under the same assumption on 13(e)" 1 as Theorem 4.3, 
^(t) := fi^(^(t); •) weakly converges to fi(t) := p(^ t)tV{t) (y)) dy as 
e I in the space C([0, T], VJlp) for every T > 0, where p = sup 0<£<1 e d N(e). 
The process (9(t),r](t)) is characterized by the three properties (i)-(iii) in 
Theorem 4.3 and p(y) is the macroscopic limit density function of the ini- 
tial configuration . 

Before giving the proofs of Theorem 4.3 and Corollary 4.4, we state a 
proposition whose proof will be postponed to the next section. 
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Proposition 4.5. For every v > 0, there exist C > and < £o < 1 
such that 

de 

— T (x) 



(4.13) 



(4.14) 



dxj 



de 99 

( X )-7JI7( Z ( X )) 



dx] dx] 



<Ce K ~\ 
< Ce d+U+1 , 



hold for everyxeM^'ie^Aciz^)), 1 < 7 < d, 1 <i< N(e) and < e < 
e . 

Proof of Theorem 4.3. Step 1. Prom the SDE (3.9) for x (e )(t), since 
V£7(x) = — VU(— x) holds from the radial symmetry of U, we have 

_l-«/2 

(4-15) ^(t) = 

i=l 

which is equivalent to e 1 ~ K / 2 A^ _1 / 2 t(;(t) in law. This shows the property (ii) 
for the limit rj(t) ofrj £ (t) noting that p = X\m e iQE d N so that lim e jo e l ~ K ' 2 N~ 1 ' 2 

p~ 1/2 - 

Step 2. We next consider the limit of 9^ e \t). Let er^ be the stopping time 
defined by (3.5) with c(e) satisfying the conditions in the theorem. Then, 
9{x^ £ \t)) is well defined for t < o-( £ \ and again from the SDE (3.9) and using 
the property (4.2) of 0(x), we have 

N 

(4.i6) de( e \t) = £- K / 2 J2v Xi e(x (E Ht))°dwi(t), t<a^\ 

1=1 

by applying Ito's formula, where V x .9 o dwi := £L=i 99 /dxj o dto^ G M{d). 
This may be further rewritten as 

d0< e >(t) = 0( £ )(t) o dm^i), t < ff ( e ), 

with an so((i)-valued martingale m^ s \t) = {m a ^'^ (t)}^ g =1 defined by 

N ft 

m^{t) = e^ 2 Y j 9{^\s))- 1 V Xi 9{^\ S ))dw i {s), t<a&. 

Therefore, setting 9&(t) :=9^(tAa^) and m^ £ \t) :=m^(tAa^) for all 
t > 0, 9^ (t) satisfies an SDE 

d9^(t) = 9^(t)odm ( - £ \t), t>0. 

This SDE written in Stratonovich's form is equivalent to 

d^ £ \t)=9 {£ \t)dm {£ \t) + \^ £ \t)d{m^ £ \rh^){t) 
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in Ito's form, where the quadratic variational process 
is defined by 



((m( £ U (e) )W)^ = £ (m Q7 ' (e) ,m 7/3 ' (e) )(t), l<a,0<d. 

7=1 

The goal is to show that 9^ (t) weakly converges as e j to the solution 
6(t) of the SDE (4.12) which is equivalent to 

d9(t) = 6(t)dm(t) + \9{t)d{m,m){t) 

in Ito's form. Indeed, once this is proved, since lim £i0 P(cj( £ ) > T) = 1 from 
the microscopic shape theorem (Theorem 3.4), 9^ £ \t) also weakly converges 
to 0(t); recall that 9^{t) was arbitrarily defined after the time when x( £ ) (t) 
goes outside of (c(z^)). To show the weak convergence of 9^(t) to 
9(t) in C([0, T], SO(d)), it suffices to prove the following two conditions for 
the driving martingale rh^ £ '(t): 



(4.17) limE 

e|0 



E[{m a ^ £ \m ab <^){t)\F s 



t(6 aa 5 l3b - rt^") 







q a + q p 

for every t > s > 0, 1 < a, (3, a, b < d, where J- s = a{w(s');s' < s} and 

d 



(4.18) 



sup sup E 

0<£<10<t<T 



dt 



{m a ^ £ \m ab ^){t) 



< oo 



for some p>2; see, for instance, [11], page 222, Theorem 5.2.1. Note that, 
from the property (iii), the quadratic variational processes of (m a ^(t)) a g are 
given by 



{m al3 ,m ab )(t) 



t(5 aa 5 l3b - 



q a + qP 

Now let us prove (4.17) and (4.18). Since we have 
^-(fh a ^ £ \rh ab ^)(t) 



dx- 



i=l 7=1 

for t < cr& (and = for t > c^), (4.13) in Proposition 4.5 implies 
d 



(4.19) 



dt 



{m a ^ £ \m ab ^){t) 



<Ce~ K N(e K - 1 ) 2 <C 
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which is bounded in e and t. This shows (4.18). To prove (4.17), define 

m( £ Xt) by 



(4.20) m( £ )(i) = e" K/2 V / 9(z(^ £ \s))y 1 V Xi 9(z(^ £ \ S )))dw i ( i 



We have replaced xS e \s) with z{pc^ e \s)) in Then, since condition 1 

on implies R(z^) < Roe" 1 for some Rq > 0, Lemma 2.4 shows xS E \s) G 

7WS (£) (c(e)) C M&(R(z^)c(e)) C ^(Hoe* -1 A c(z^)) c .M^V -1 A 
c(z( £ ))) for s < <j( e ) by taking v' G (0, i/) and for small e > 0. Therefore, by 
(4.14) in Proposition 4.5 and noting that 0(x) = 0(z(x)), we have 

(4.21) |(m a/3 '( £ ) - m Q/3 '^)(i)| < Ce- K N(e d+v ' +1 ft < C'e 2u 't. 
Using (4.19) and (4.21), we get 

\(m a/3 '( £ \rh ab ^)(t) - (m a ^ £ \fh ab '^){t)\ 



< yj (m«A( £ ) - fh a P>( e ))(t) • (m a6 '( £ ))(t) 

+ \J {m ab ^ £ ) - m afe -( £ ))(t) • (rh a P^)(t) 

< C"e u 't, 

which tends to as s [ 0; note that a similar estimate to (4.19) can be 
shown for (fh a/3 '^ £ \m ab '^)(t). We accordingly see that it suffices to show 
condition (4.17) for fhS e > instead of rhS £ \ 

For a configuration z = (fg -(z^) £ with some (6, fj) G SO(d) x 

since 9(z) = 9 and Q(z) = Q{z^)9~ l , Proposition 4.2 shows that 



09 



( n Jz&)) = 0{ Proj o(Q(z( £ )))}" 1 Proj {(e~\) ® zf} 



However, we see that 

(4.22) lime K Q(z^) =limf / y a y^^ £ \z^; dy)) =Q 

ej.0 \JM.d / a,/3=l 

(4.23) {Proj 0(0)}"^ 



for Y = (y a(S ) a(3 eso{d), and 

proj {(*-%) gi = i(^zf' (e) - ^^r ,(£) )^=i 5 

where 9 = {9 a ^) a p and = (zf ) a - Therefore, from Lemma 5.1(h) stated 
below [note that {Proj o(Q(z^))}~ 1 = e K ^ £ \z^)~ 1 ] and recalling e\zf'\ < 
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R, we have 



(4.24) lim 

e|0 



*»6-^<u>- (=<*)) J Q^' ie) -e^?' {£) \ d 



a,/3=l 



uniformly in (9,fj). Hence, 



N d 



(^ (e) (t)2?' (e) - n p M{t)z? {e) ) 



x {6^ £ \t)z\' {£) - 9^ £ \t)Zi' [E) ) + o(l) 

6 aa §l 3b _ S ab 5 f3a 



q a + q' 3 



+ o(l), i < cr (e) , 



as e I for every 1 < a, /3, a, b < d. The error terms o(l), which come from the 
errors in (4.24), tend to uniformly in t and u [an element of the probability 
space on which frS £ \t) are defined]. This proves (4.17) for fh( £ \t). 

Step 3. Finally, to show property (i), compute the quadratic variational 
processes of m afi ^ £ \t) and rp'®(t) from (4.15), (4.20) and (4.24): 

^-(m a ^ £ \r]^)(t) 

= 4tY- — -* — + o(l), t<a( £ \ 

However, since is centered, the sum vanishes and this proves property (i). 
□ 



Proof of Corollary 4.4. Theorem 4.3 combined with Theorem 3.4 
shows that (f,fi^ £ \t)) weakly converges to (/, /x(t)) as e j in the space 
C7([0,T],K) for every / G C b (R d ), where (/,//) := f Rd f(y)fi(dy). Therefore, 
to conclude the corollary, it suffices to show the tightness of the family 
of laws of {/i (e )(i); < e < 1} on the space C([0,T],Wlp). But this can be 
deduced from: 

1. For each 5 > 0, there exists a compact set K, in 9Jtp such that P{n^ (t) G K, 
for every te[0,T])> 1-6. 

2. For every / E C b (R d ), {{f,^ £ \t)); < e < 1} is tight on the space C([0,T],M). 

Condition 2 follows from what we mentioned above. Condition 1 is also easy, 
since the support of ^ £ \t) is in the ball with center r/ £ )(i) and radius R 
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as long as t < o~^ e ' and Tj^ e '(t) ==> Tj(t), which is the time changed Brownian 
motion, as e j 0. □ 

The solution 9(t) of the SDE (4.12) is called the left Brownian motion on 
SO(d); see [15, 16]. A coordinate satisfying the relation like (4.2) was used 
by Katzenberger [9] to make a cancellation for diverging terms as we have 
seen in deriving (4.16). 

Example 4.1. (i) Let D be a bounded domain in R d (d = 2 or 3) hav- 
ing a smooth boundary dD and let z^ be the infinitesimally rigid crystal 
constructed from e _1 D n aA on the d-dimensional triangular lattice as in 
Example 2.3. The configuration is the microscopic crystal consisting of 
atoms arranged in an equal distance a and D is the corresponding macro- 
scopic body. The macroscopic density function of z^ is given by 



p(y) 



My) 

- d |detA|' 



where A = {e\&i ■ ■ ■ e^) G M(d) is the matrix consisting of d column vectors 
{e a } a used for the definition of the triangular lattice. In this sense, D is the 
high density region and the outside of D is the empty region. 

(ii) In higher dimensions, one can construct z^ based on the idea ex- 
plained in Example 2.1. 

5. Proof of Proposition 4.5. This section gives the proof of Proposition 
4.5. Consider an operator <I> on the space so(d) defined by $>X = Proj(<5A) 
for X £ so(d). Then, as we have seen in (4.23), <& is invertible and the 
operator norm of $ _1 can be dominated by 

(5.1) H* -1 !!^^^ max 



Recall that Q(x) = Q( £ )(x) is determined from z^ and set Q( e \~x) = e K Q^{x) 
Then, we have Q = lim e j_oO^( z ^)> see (4.22). We introduce another oper- 
ator $( £ )(x) on so(d) by 

$ (e) (x)A = Proj (0 (£) (x)fl(x)X), X e so(d). 

In the following, we denote M^(e v - 1 Ac(z< £ ))) simply by M^^' 1 ), since 
we only use the bound ||h(x)||oo < e u ~ l for x. 

Lemma 5.1. (i) For every v > 0, there exists Eq > such that 
sup sup ||$^(x) _1 || < oo. 

°<e<eo x6M W( B ^-l) 
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(ii) For every v > 0, 



lim sup IIS^x) -1 -^ -1 !! =0. 



Proof. We first assume that x e Moo *) satisfies z(x) = z^ s \ Then, 
we have 



(5.2) 



|Q (£) (x)-Q|<r(e), 



with some r(e) — ► as e J 0. Indeed, the left-hand side of (5.2) is dominated 
by 

|Q( £ )(x) - Q( £ )(z( £ ))| + \Q^(z^) - Q\, 
and the first term is further bounded as 



N 



\e K q a P(x)-e K q al3 (z&)\<e K Y^ 



a,(e)i|J Me) 



i=l 



Xi — Z; 



< e K ■ N ■ Re~"e v ~" < Ce u , 

while the second term tends to as e J 0. 

Denoting by ^ = <3?( £ )(x) — $, <3>( £ )(x) _1 can be expressed as 

oo 

^(x)- 1 = (j+ $- i ^)- i $- 1 = ^(-s- 1 *)*^- 1 . 

k=0 

Since 0(x) = 6(z^) = I, we have = ±{<3 (e) ( x ) - Q}X + ±X* {Q& (x) 
Q} and therefore \\^f\\ <r(e) from (5.2). Accordingly, (5.1) implies 



5^(-$- 1 *) fc 

fc=0 



< 



1 



1 - Cr(e) 



<2, 



so ||$( e )(x) _1 || < 2C for sufficiently small e > such that r(e) < 1/(2C). 
By acting rotation and translation, similar estimate can be derived for 
||^>^> (x) _1 1| for every xeM^e"" 1 ) [without assuming z(x) = z^] and 
this concludes the proof of (i) . The second assertion (ii) follows from 



fc=l 



< 



Cr{e) 



1 - Cr(e) 



□ 
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PROOF of (4.13) in Proposition 4.5. By Proposition 4.2, 
|^(x) = e^(x)$( £ )(x)- 1 Proj {^(x)- 1 ^) ® 
and therefore (4.13) follows from Lemma 5.1 (i) noting \zf \ < Re -1 . □ 



Lemma 5.2. For ewery 1 <j,j' < d, 



sup sup sup e 

0<K £ o xeA1 g {£ ,-l ) Ki,i'<N( E ) 



-2k+2 



d 2 e 



Ox] Ox], 



7(X) 



< oo. 



Proof. The identity (4.9) may be rewritten as 



06 



Q^(x)^(x) + (4 e) ®e 7 )fl(x),y) 



Hence, taking the derivative of this identity in xj, , we have 

a 2 t 



Q (£) (x) 



Ox] Ox- 



7 



-(x) + (zf 9 ey)|^(x) + (zf J ®g^(x),y)= 0. 

"^t oxL / 



This implies 



Ox] Ox], 



^(x)$( £ )(x)~ 1 



, M . 06 , , .(e) . 86 . . 
x Proj ^ (4 £j ® eyj^-y (x) + 8) e 7 )^-y (x) 



The conclusion follows from Lemma 5 . 1 (i) and (4.13) by noting that \z#\, 
|4 e) |<ite _1 . □ 

PROOF of (4.14) in Proposition 4.5. Applying the mean value theo- 
rem, we have, from Lemma 5.2, 



06 



(x) ™(«(x)) 



N d 

<££ 

i'=iy=i 



Ox] Ox], 



t(x*) 



1*2 - (z(x))7, 



<CiV-e 



where x* is a certain point on the segment connecting x and z(x). This 
shows (4.14). □ 
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6. Concluding remarks. 

6.1. We have studied in Section 4 the case where the macroscopic body 
is d-dimensional, but one can consider thin bodies and derive their motion as 
well. Assume that, for an n-dimensional Riemannian manifold M in M. d ,n < 
d, a sequence of infinitesimally rigid crystals z& = ){Li is given and it 
has a macroscopic limit density function pm{u) on M in the sense that 

N 

(6.1) lime n ^<5 u){dy) = p M {y)dy M , 

i=l 

where dyu denotes the volume element of M. Then, comparing (6.1) with 
the condition 2 for in Section 4.2, we see that p(y) = for the sequence 
z( e \ which we are considering here, under the scaling (4.10), and therefore 
p = q a = in Theorem 4.3. This means that a different time scaling is re- 
quired for x(i) to have a nontrivial macroscopic limit. Indeed, one can show 
that the right scaling is x^ £ \t) := x(e~( n+2 )t), and under this scaling, ran- 
dom motion of the body (M,pm) is obtained in the limit. Note that the 
affine hull of z[ £ > should be d-dimensional to be infinitesimally rigid (see [2], 
page 174) and therefore, even for obtaining an n-dimensional macroscopic 
body in the limit, ci-dimensional configurations should be considered mi- 
croscopically. High polymers or membranes studied in physical chemistry 
usually have the above structures with n = 1 or 2 in IR 3 . 

6.2. If the Hamiltonian is suitably modified, the notion of rigidity for 
the microscopic configurations may change. Let V E Cg(M) be a symmet- 
ric function having a deep well at satisfying V"(0) > and consider the 
modified Hamiltonian of H(x) by adding a three-body interaction term: 



iT(x)=if(x) + £V 



%i ~\~ Xj 



Xk 



{diam{xi,a;j,a;j.}<ao} ; 



where the sum is taken for i, j, k different and ao is a constant smaller than 
4a. We assume b <2a for the potential U in H(x). Then, one-dimensional 
straight chains z = {zi)f =l in M. d :zi — zq = i{z\ — zq), 1 < i < N, arranged in 
an equal distance a (i.e., \z± — zq\ = a) are local minima of H. The proper 
time change is x^ e \t) := x(e~ 3 t), which is the same as taking n = 1 in Sec- 
tion 6.1. The particles' number iV should behave as p = lining exists. 
Introducing an energy different from ours, Kotani and Sunada [10] charac- 
terize the equilibrium configurations of crystals. 
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